#4.2.1.

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")

summary(cervejas)

cervas <- table(cervejas)
cervas
plot(cervas)

cervass <- as.numeric(cervas)

dotchart(cervass)

#4.2.2. O grfico de barras.

#4.3.1

caixeta <- read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=dados:caixeta.csv", sep = ",", header = T)

raio <- caixeta$cap/(2*pi)
DAP <- raio*2

hist(x = DAP, breaks = 40, xlim = c(0,350))


#4.3.2

hist(x = caixeta$h[caixeta$local=="chauas"],  breaks = 50, xlim = c(0,200))


hist(x = caixeta$h[caixeta$local=="jureia"],  breaks = 50, xlim = c(0,200))


hist(x = caixeta$h[caixeta$local=="retiro"],  breaks = 20, xlim = c(0,200))

#grafico unico para vrios locais
histogram(h~local, data = caixeta)

#4.3.3

#Sim, a distribuio de tamanho  diferente entre os locais.

#4.4.1

e.grandis <- read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=dados:egrandis.csv", sep = ";", header = T)

hist(x = e.grandis$dap[e.grandis$regiao], breaks = 4, xlim = c(6,12))

hist(x = e.grandis$dap[e.grandis$rotacao],  breaks = 4, xlim = c(6,10))

#4.4.2

hist(x = e.grandis$dap, breaks = 20, xlim = c(0,40))
qqnorm(e.grandis$ht); qqline(e.grandis$ht)

#4.5.1

caixetaa <- cbind(caixeta,DAP)
tabebuia <- caixetaa[caixeta2$especie=="Tabebuia cassinoides",c("DAP","cap","h","arvore","local","especie")]
areabasal <- pi*((tabebuia$DAP/2)^2)
tabebuiaaa <- cbind(tabebuia,areabasal)
tabebuiaaa <- aggregate(x=tabebuiaa$areabasal,by=list(arvore=tabebuiaa$arvore,local=tabebuiaa$local,h=tabebuiaa$h),sum)
colnames(tabebuiaaa) <- c("arvore","local","altura","abt")
daptotal <- 2*(sqrt(tabebuiaaa$abt/pi))
tabfinal <- cbind(tabebuiaaa,daptotal)

par(mfrow=c(1,3))
with(subset (tabfinal, local == "chauas"), plot (daptotal~altura, xlab = "DAP", ylab = "Altura", main = "Chauas"))
with(subset (tabfinal, local == "retiro"), plot (daptotal~altura, xlab = "DAP", ylab = "Altura", main = "Retiro"))
with(subset (tabfinal, local == "jureia"), plot (daptotal~altura, xlab = "DAP", ylab = "Altura", main = "Jureia"))

#4.5.2
par(mfrow=c(1,3))
with(subset(tabfinal, local=="chauas"), scatter.smooth(daptotal~altura, col="blue", xlab="DAP", ylab="Altura", main = "Chauas"))
with(subset(tabfinal, local=="retiro"), scatter.smooth(daptotal~altura, col="blue", xlab="DAP", ylab="Altura", main = "Retiro"))
with(subset(tabfinal, local=="jureia"), scatter.smooth(daptotal~altura, col="blue", xlab="DAP", ylab="Altura", main = "Jureia"))
 
#4.5.3
library(lattice)
xyplot(daptotal~altura|local, data=tabfinal, xlab="Altura", ylab="DAP")